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Abstract 

The Bose-Einstein condensation (BEC) critical temperature in a relativistic ideal 
Bose gas of identical bosons, with and without the antibosons expected to be pair- 
produced abundantly at sufficiently hot temperatures, is exactly calculated for all boson 
number-densities, all boson point rest masses, and all temperatures. The Helmholtz 
free energy at the critical BEC temperature is lower with antibosons, thus implying 
that omitting antibosons always leads to the computation of a metastable state. 

PACS # 03.75.Kk; 05.30.Jp; 47.75.+f 

Since its theoretical prediction by Einstein in 1925 based on the work in 1924 by Bose on 
photons, and after languishing for many decades as a mere academic exercise in textbooks, 
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BEC has been observed in the laboratory in laser-cooled, magnetically-trapped ultra-cold 
bosonic atomic clouds of gfRb [1J, ^Li [2J, ffNa [3], \R @], ffRb 0, pie [BJ, [7], |fCs 
[5] (a), 7Q 4 Yb [H](b), and ^Cr [9]. Previously, BEC in a gas of excitons had also been reported 
[ID] . More recently, BEC has been seen as well in fermionic atomic gases of fgK [IT] and 3LI 
[12] as a result of some of the fermions presumably Cooper-pairing [13] into bosons. The 
role of hole Cooper pairs accounted for or not, along with the usual particle Cooper pairs, 
has been explored [H]-[I7] with striking implications for many-electron superconducting 
states within a BEC scenario. Their effects in neutral-fermion many-particle systems are 
yet to be investigated. It is the analogy of hole pairs with antibosons which has piqued our 
interest in the problem to be dealt with in this paper. Although it is hard to imagine a 
many-boson system whose constituent bosons do not disintegrate into their components at 
the high temperatures where antiboson production becomes substantial, it may shed light 
on the dynamics of many-fermion systems where both particle- and hole-Cooper-pairing can 
occur at all temperatures. Moreover, recent specific-heat measurements [18] in TlCuCl^ 
suggest a magnon BEC with a relativistic dispersion relation. 

In early papers [19]- [21] on the relativistic ideal boson gas (RIBG) explicit BEC criti- 
cal transition temperature T c -formulae were derived for both the nonrelativistic and ultra- 
relativistic limits and specific-heat anomalies at T c were studied. In addition, Refs. [20j[21j 
considered all space dimensions d > and delved into the relation between d and various 
critical exponents. Antiboson production, however, was not accounted for. The first papers 
to include both bosons and anti-bosons appear to be Refs. [22] [23] where high-temperature 
expansions for the various thermodynamic functions (pressure, particle-number density, en- 
tropy, specific heats, etc.) were derived. Extensive numerical work in d dimensions that 
does not rely such high-temperature expansions was reported in Refs. [21] [2S] • In an elegant 
treatment [26] with inverse Mellin transforms the specific heat anomaly of the RIBG at its 
BEC T c was found to be washed out when pair-production was included. The relationship 
between the BEC of the RIBG and spontaneous-symmetry breaking was explored in Refs. 
[23j[27J; see also the rather complete Ref. [28J, esp. Sec. 2.4. The so-called BCS[29]-to-Bose 
crossover scenario (see Ref. [17] and refs. therein), and even the pseudogap concept [30] of 
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superconductors, first seemed to have appeared in quark physics in Ref . [31] ; for a review see 
Ref . [32] • More recently, two "crossovers" have been identified [33] in an interacting fermion 
gas where pairing into bosonic Cooper pairs [13] can occur to form a relativistic superfluid 
as an example of the BCS-Bose crossover followed by a Bose-to-RIBG /BEC-crossover where 
both anti-bosons as well as bosons dominate the thermodynamics. A fully relativistic detailed 
study [31] of these crossovers at zero temperature has also appeared. 

In this Letter we exhibit, as a function of boson number density, exact BEC transition 
temperatures for the RIBG gas of identical bosons with and without antibosons in 3D. The 
system with both kinds of bosons always has the higher T c , i.e., is the system with the first 
BEC singularity that appears as it is cooled. This suggests that the Helmholtz free energy 
might be lower and thus correspond thermodynamically to the stable system as opposed 
to a metastable system for the lower-T c system. It is then calculated and indeed found to 
be lower at all densities for the complete problem with both bosons and antibosons, when 
compared to the problem without antibosons. This implies that the omission of antibosons 
will not lead to stable states. 

The number of bosons N of mass m that make up an ideal boson gas in d dimensions 
(without antibosons) is iV = Xlk n k = ^2[e^ Ek '~^ T '^ — where (3 = 1/kgT, ks is the 

~" k 

Boltzmann constant and /~t(T) is the boson chemical potential. Here, the total energy of 
each boson is 

| E k | = Vc 2 h 2 k 2 + m 2 c 4 (1) 

= mc 2 + h 2 k 2 /2m + 0(k 4 ) if chk < mc 2 NR (2) 

= chk[l + ^(mc/hk) 2 + 0(k~ i )} iichk^>mc 2 UR (3) 

where k is the boson wavenumber, m its rest mass and c is the speed of light. The two limits 
refer to the nonrelativistic (NR) and ultrarelativistic (UR) extremes. For a cubic box of 
side length L in the continuous limit the sums over the <i-dimensional wavevector k become 
integrals as ~ * (L/2n) d J d d k. At the BEC critical transition temperature T c , /x(T c ) = mc 2 

k 
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and the boson number density can be expressed as 

n = — = — / d d k - (4) 

~ L« (2vr) d 7 exp[f3 c (\E k \-mci)}-l [) 

where (3 C = l/k B T c . In the nonrelativistic (NR) extreme ([2]) inserted into (j4]) leaves 
n = (2ir)~ d J d d k [exp [/3 C (h 2 k 2 /2m)] - l]" 1 . Putting d d k = [2ix d l 2 /T(d/2)}k d - 1 dk when 
integrating over terms independent of angles gives an expression for n as function of T c in 
terms of the Bose function [35] g a (z) of z = exp(fi/k B T c ) which for z = 1 diverges, namely 
g a {l) — > oo when a < 1, but becomes the Riemann Zeta function £(er) < oo when a > 1. 
Here r(cr) is the gamma function. Solving for the critical temperature then gives 

k B T c NR ~ B = ^ [n/C(d/2)f d (5) 

where the superscript NR-B stands for the nonrelativistic limit with bosons (B) but no 
antibosons (B). In 3D this reduces to the familiar textbook result k B T^ R ~ B ~ 3.31h 2 n 2 ^ 3 /m 
since C(3/2) — 2.612. In the ultrarelativistic (UR) extreme the leading term of ([3]) inserted 
into (jlj) leads to T c = for all d < 1 since then <?<i(l) diverges. However, for d > 1 

~Ti d c d 2 d - 1 it d l 2 T{d/2y 1/d 



k B T^ R - B 



n 



l/d 



(6) 



r(d)C(d) 

which in 3D becomes k B T^ R ~ B = hc7i 2/3 [n/((3)] 1/3 ~ 2.017W/ 3 as C(3) ~ 1.20206. In 2D 
T^ R ~ B ^ unlike the common instance with quadratic dispersion where T c vanishes because 
g x (l) diverges; specifically k B T^ R ~ B = hc[2nn/C(2)} 1 / 2 ~ 1.954/icn 1 / 2 since C(2) = vr 2 /6. 

At sufficiently high temperatures such that k B T 3> mc 2 boson- antiboson pair-production 
occurs abundantly; this has been stressed by Huang [36] . The total energy E k of each parti- 
cle always satisfies E\ = c 2 h 2 k 2 + m 2 c 4: so that E k — ±\E k \ where | E k | is given by (JTJ and 
with the + sign referring to bosons and the — sign to antibosons. Instead of iV = J^ k ^k 
the complete number equation is now [22] 

N _ n = £(n k -n£) (7) 

k 

1 1 



E 



k 



exp [/3(| E k | -//)]- 1 exp [/3(| £ fe | + //)]- 1 
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where n^, (n^) is the average number of bosons (antibosons) in the state of energy ± | E\. \, 
respectively, at a given temperature T and N (N) is their respective total number at that 
temperature. Since n^, > for all k and Eq = mc 2 , the chemical potential must be 
bounded by —mc 2 < fi < mc 2 . Instead of N constant one must now impose the constancy of 
N — N to extract the correct BEC critical temperature, say, T B B referring to both bosons 
(B) and antibosons (B). Since fi(T BB ) = mc 2 (J7|) becomes 

oo 

/ -jL-r ~tCt\ i t (i 2ir d/2 f „ . , sinh( (3 c mc 2 ) 
n = (N- N)/L = — -r dkk^ 1 v — J- . (8) 



r(rf/2)(2vr) d J CQsh (j3 c y/ c 2%2 k 2 + m2c 4^ _ cos h(/? c mc 2 ) 

This is an exact expression for the BEC T c of an ideal Bose gas at any temperature as it 
includes both bosons and antibosons; it is consistent with Eq. (13) of Ref. [24J. 

At low enough temperatures such that ksT c <C mc 2 antibosons can be neglected and 
(JSJ) simplifies to, say, T C 7VB " BB which is precisely (JSJ) as expected. In the opposite extreme 
ksT c ^> mc 2 ([8]) leads to the limiting expression, say, 



k B T™- BB 



n d c d ~ 2 T(d/2)(2>n) d ll/(d ' 1] 



n 1 '^ (9) 



_Am7c d / 2 T(d)C{d - 1)_ 

that sharply differs from (jHJ). In 3D this becomes ksT^ R ~ BB = (Sh^c/my^n 1 / 2 , a result ap- 
parently first reported in Ref. [22]. This novel relation has suggested [18] itself experimentally 
as a magnon BEC in specific- heat measurements in TlCuCl^. 

As functions of the dimensionless boson number-density h 3 n/m 3 c 3 Fig. 1 displays the 
behavior of the exact T BB (in units of mc 2 //cs) numerically extracted from (jSJ) for d = 3 
(thick full curve labeled "exact BB") compared with the nonrelativistic limit T^ R ~ BB from 
(jSJ) (dashed line labeled "NR-BB" ) and with the ultrarelativistic T^ R ' BB just stated (full thin 
line labeled "UR-BB"). Fig. 2 shows how, at sufficiently high densities n and/or sufficiently 
small boson rest mass m, the exact T BB (again in units of mc 2 /^B, full curve labeled "exact 
BB" ) is clearly the first BEC singularity encountered as the many-boson system is cooled, 
compared with the "later" BEC singularity in the system without antibosons at T B (dashed 
curve labeled "exact B") extracted numerically from pi) . It is then tempting to speculate 
that the boson gas with both kinds of bosons will be the more stable, i.e., have a lower 



Helmholtz free energy at all critical temperatures, at any fixed h 3 n/m 3 c 3 . This will now be 
shown to be the case indeed. 

The exact Helmholtz free energy per unit volume V = L 3 for the boson-antiboson 3D 
mixture, when T = T c = l/ksPc and \x = mc 2 , is 

POO 

F exactBB (T c , V)/V = nmc 2 + {k B T c /2n 2 ) / dkk 2 {\n[l - exp(&[mc 2 - | E k |])] 

Jo 

+ln[l - exp(-/? c [mc 2 + | E k |])]} (10) 

where | | is given by ([1]) and n = (N — N) /V, and T c is extracted numerically from (jSJ) 
for each value of n. If ksPc << mc 2 antibosons can be neglected entirely and F exactS (T c , V) 
is just OH]) but without the second log term and for which T c is now extracted numerically 
from (j3J) instead of from (jSj) for each value of n. Figure 3 is the difference between the 
two free energies and more clearly shows why F exactBB (T c ,V) is always lower. The inset 
figure shows the behavior of these two free energies, each evaluated at their appropriate T c 
value; by inspection, both curves correspond as they should to positive pressures P since 
P = —(dF/dV)N,T- Figure 3 proves the speculation advanced. 

In summary, based on exact numerical calculations of BEC T c s in a 3D RIBG with and 
without the antibosons expected to be pair-produced at higher and higher temperatures, 
the highest critical T c is that associated with the RIBG system with both kinds of bosons 
taken into account. At lower n and/or larger m the higher T c merges smoothly from above 
onto the lower T c of the RIBG system with antibosons neglected. Comparing the associated 
Helmholtz free energies shows that the RIBG with both kinds of bosons has lower values 
for all n and m and thus substantiates the initial suspicion that the RIBG system with no 
antibosons is metastable with respect to the one with both kinds of bosons. 
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Figure 1. BEC T c s (in units mc 2 /kg) as function of boson number- density n 
expressed in dimensionless form as h 3 n/m 3 c 3 . Thick curve labeled "exact BB" is 
exact numerical result of ([H]) that corresponds to BEC T c in a RIBG with both 
bosons B and anti-bosons B. Thin full straight line labeled "UR--B-B" is the 
ultrarelativistic limit (Q for d = 3 with both kinds of bosons. Dashed straight 
line labeled "NR-.B.B" is its corresponding nonrelativistic limit ([5]) for d = 3 and 
tends asymptotically to the "exact BB" curve at smaller n. 



Figure 2. Same as Fig. 1 comparing exact BB RIBG T c extracted from ([8]) 
against that of exact B from (Tj0) . 
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Figure 3. Difference between exact Helmholtz free energy density F/V (in 
units of nmc 2 which is the total rest-mass energy density) (full curve labeled 
"exact.B.B" ) and that without anti-bosons (dashed curve labeled "exact-B" ) using 
same horizontal axes as in Figs. 1 and 2. The inset figure displays the two free 
energies. 

under Contract # DE-AC 52-06NA25396. 
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